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Abstract:  Breakdown  of  a  slender  longitudinal  vortex  caused  by  a  normal  shock  is  studied  with 
a  numerical  solution  of  the  Euler  and  the  Navier-Stokes  equations  for  time- dependent,  three- 
dimensional,  supersonic  flow  at  a  free  stream  Mach  number  of  1.6  and  different  vortex  strengths. 
When  breakdown  occurs,  a  free  stagnation  point  is  formed  downstream  from  the  shock,  followed  by  a 
region  of  reversed  flow  with  a  bubble-like  flow  structure.  The  burst  part  of  the  vortex  grows  in  the  axial 
and  radial  directions  until  the  bubble  reaches  a  stable  position.  The  shock  remains  normal  near  the 
axis  of  the  vortex  but  is  curved  further  away  from  it.  The  flow  downstream  from  the  shock  is  slightly 
oscillating.  The  numerical  results  clearly  reveal  the  time- dependent  flow  structure  in  the  axial  and  the 
radial  direction.  The  results  compare  well  with  recent  experimental  findings. 

1  Introduction 

In  continuation  of  previous  studies  [1,  2,  3,  4]  the  interaction  of  a  longitudinal  Burgers  vortex  with  a 
wake-like  axial  velocity  profile  and  a  normal  shock  is  investigated  for  viscous  supersonic  flow  at  a  Mach 
number  of  Ma=1.6.  Such  an  interaction  is  encountered,  for  example,  on  st rake- wing  configurations  at 
moderate  and  high  angles  of  attack  and  compressors  operating  near  their  stability  limits  [5],  where 
the  tip  vortex  of  the  blades  interacts  with  a  normal  shock  during  the  passage  through  the  compressor 
cascade.  Although  the  breakdown  process  is  mainly  dominated  by  pressure  forces,  the  flow  inside 
the  burst  part  of  the  vortex  is  also  strongly  influenced  by  viscous  forces.  A  comparison  of  results  of 
previous  computations  for  inviscid  flows  with  those  for  viscous  flows  presented  here  clearly  exhibits 
that  the  vortex  structures  encountered  in  viscous  flows  inside  the  bubble  are  smaller,  and  that  the 
bubble  itself  is  smaller  in  its  radial  and  axial  extent. 


2  Governing  equations,  boundary  conditions,  and  computa¬ 
tional  domain 

The  numerical  solution  used  in  the  present  analysis  is  based  on  the  Euler  and  the  Navier-Stokes 
equations  for  time-dependent,  three-dimensional,  compressible  flow.  The  convective  terms  are  dis¬ 
cretized  with  an  explicit  finite- volume  scheme  based  on  the  AUSM+  approach  by  Liou  [6].  The 
dissipative  terms,  including  the  Stokes  stresses  and  the  Fourier  heat  flux  vector,  are  discretized 
using  an  adapted  cell- vertex  scheme,  and  the  time-derivatives  are  approximated  via  a  5-step  Runge- 
Kutta  scheme.  Details  of  the  computational  method  are  presented  in  [2].  The  computational  domain 
is  a  rectangular  box  with  99x99x199  grid  points  in  the  x-y -,  and  ^-direction  (Fig.  1).  Near  the  axis 
of  the  vortex  the  grid  is  clustered  so  that  the  vortex  core  (r  <  1)  is  resolved  with  approx.  21  grid 
points  in  the  radial  direction.  The  inflow  and  initial  conditions  are  prescribed  with  a  slender  Burgers 
vortex  in  dimensionless  form 

vv(.r)  = ■  eL^~  ;  vz(r)=vZoo-(l-S-e~^-r2)  ;  vr{r)  =  0  (1) 
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In  the  above  equation  v^(r)  is  the  circumferential,  vz(r)  the  axial,  and  vr(r)  the  initial  radial  velocity 
component.  Figure  2  shows  the  vortex  parameters  as  a  function  of  the  core  radius  at  a  free  stream 
Mach  number  of  Maoo  =  1.6.  The  quantity  T0  =  To/(ao  *ro)  is  the  circulation  of  the  vortex  scaled  with 
the  speed  of  sound  at  stagnation  conditions,  a g  and  the  radius  of  the  vortex  core  at  r  =  r*/rg  =  l. 
The  quantity  (5=0.1  determines  a  wake-like  axial  velocity  profile,  and  r  =  ±y^l/(2 fiw)  describes  the 
position  of  the  point  of  inflection  for  fi  =  1.  The  initial  radial  pressure  distribution  is  obtained  by 
solving  the  radial  momentum  equation  for  slender  vortices;  the  temperature  is  computed  by  assuming 
constant  stagnation  enthalpy  for  the  radial  direction,  and  the  equation  of  state  yields  the  density.  For 
the  subsonic  part  of  the  outflow  plane  (Fig.  1,  plane  6)  a  non-reflecting  boundary  condition  based 
on  the  characteristic  approach  by  Poinsot  &  Lele  [7]  for  viscous  flows  is  chosen.  The  shock  is 
prescribed  by  the  Rankine-Hugoniot  relations  in  the  x-y  plane  at  z= 0.  The  remaining  horizontal 
and  vertical  boundaries  of  the  computational  domain  (Fig.  1,  plane  1  to  4)  are  represented  by  stream 
surfaces,  so  that  the  kinematic  flow  condition  holds.  The  computation  was  initialized  by  prescribing 
the  shock  in  the  middle  of  the  computational  domain.  The  main  flow  is  in  the  ^-direction.  Upstream 
of  the  shock  the  slender  vortex  is  embedded  in  the  main  flow  whereas  downstream  a  uniform  flow  is 
specified. 


3  Results 

It  was  shown  in  [2],  that  normal-shock  vortex-interaction  does  not  necessarily  lead  to  vortex  break¬ 
down.  In  fact,  when  the  vortex  strengths,  e.g.  its  circulation,  is  small,  the  vortex  can  pass  the  shock 
without  generally  changing  its  shape.  Two  major  statements  could  be  made:  First,  by  passing  the 
shock  the  axial  velocity  profile  of  the  vortex,  which  was  uniform  upstream  of  the  shock,  becomes  wake 
like.  The  deceleration  of  flow  has  its  maximum  on  the  vortex  axis.  Second,  the  vortex  circumferential 
velocity  component  remained  nearly  unchanged  when  passing  the  shock.  The  first  can  be  explained 
by  the  pressure  minimum  on  the  vortex  axis,  the  second  finding  agrees  with  the  shock  relations.  These 
observations  were  also  made  by  Delery  [10]  in  his  experiments. 

Increase  of  the  vortex  strengths,  say  the  circumferential  velocity,  or  the  shock  strength  through  the 
free  stream  Mach  number,  or  the  vortex  deficit  of  the  axial  flow  the  interaction  with  the  normal 
shock  leads  to  vortex  breakdown.  This  can  be  seen  in  Figure  3.  The  left  column  shows  the  Mach 
number  distribution  and  the  right  column  shows  the  streamlines  in  the  symmetry  x-z-plane  of  the 
computational  domain  for  various  vortex  strengths  of  r0=3.5,  4.0,  4.5.  The  instantaneous  pictures  of 
the  Euler  computations  show  that  at  vortex  breakdown  a  free  stagnation  point  on  the  vortex  axis 
is  formed.  The  shock  is  curved  forward  and  a  bubble-like  flow  regime  grows  in  the  lateral  and  axial 
directions  until  a  nearly  stable,  slightly  oscillating  shock  topology  is  reached.  Inside  the  breakdown 
bubble  counter-clockwise  rotating  ring-like  vortical  structures  are  generated,  so  that  the  stagnation 
point  on  the  vortex  axis  is  stabilized  by  the  backflow  caused  by  the  ring-like  vortical  structures 
perpendicular  to  the  vortex  axis.  By  either  increasing  the  vortex  strength,  the  free  stream  Mach 
number  or  the  vortex  deflection  parameter  S  the  interaction  of  the  comprehensive  effects  lead  to 
strong  vortex  breakdown. 

Figure  4  shows  the  breakdown  caused  by  the  normal-shock  vortex-interaction  at  Maoo  =  1.6,  Fq=5.0, 
for  a  time  level  £=200  for  inviscid  (left  column)  and  viscous  flow  computation  (Re= 15000,  right 
column) .  The  upper  pictures  the  three-dimensional  vorticity  distribution  at  strong  vortex  breakdown 
can  be  seen.  In  these  pictures  the  Mach  number  distribution  of  the  x-z  symmetry  plane  of  the 
computational  domain  is  projected  to  the  boundary  surface  of  the  computational  domain.  It  can  be 
seen  that  the  inviscid  flow  computation  (left  column)  yields  vortical  structures  larger  than  those  of  the 
viscous  computation.  The  Mach  number  and  pressure  contours  confirm  a  multiple  shock  structure 
(Fig.  4  b).  The  deformed  shock  resembles  that  of  a  bow  shock  observed  at  supersonic  flow  around 
blunt  bodies.  Further  away  from  the  axis  the  shock  angle  decreases  until  a  shock-shock  interaction 
is  encountered.  The  calculation  for  inviscid  flow  yields  two  adjacent  counterclockwise  rotating  vortex 
rings  in  the  bubble.  For  the  viscous  computations  multiple  vortex  structures  can  be  observed  inside 
the  bubble.  An  analysis  of  the  time-development  of  the  flow  shows  a  periodic  shedding  of  the  vortex 
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rings.  A  detailed  investigation  of  the  flow  pattern  evidences  that  the  two  vortices  are  separated  by  a 
shock  which  is  caused  by  a  supersonic  jet  like  flow  in  the  upstream  direction  along  the  axis  (Figure  6). 
This  shock  is  perpendicular  to  the  axis  of  symmetry.  The  streamlines  in  the  recirculation  zone  show  a 
convergent-divergent  nozzle-like  flow  structure  that  enables  the  flow  to  be  supersonic  again  inside  the 
bubble.  The  acceleration  to  supersonic  flow  inside  the  bubble  was  already  noticed  by  Erlebacher 
et.  AL.  in  their  numerical  investigations  [9]. 

In  the  lower  part  of  Figure  4  the  Mach  number  distribution  in  a  plane  perpendicular  to  the  vortex 
axis  is  shown.  A  fan  like  flow  regime  with  multiple  shock  structures  can  be  observed  which  also  was 
noticed  by  Brillant  et.  al.  in  their  experiments  [12]. 

Figure  5  shows  a  time  sequence  from  £=ro/ao=60  to  £=200.  The  instantaneous  pictures  in  the 
two  columns  left  show  the  projections  of  the  streamlines  and  the  vortex  lines  dz/dx  =  ujz/ojx  in  the 
horizontal  center  plane  y— const,  (see  Fig.  1)  of  the  domain  of  integration.  The  third  column  represents 
a  sequence  of  A2 -surfaces,  that  to  a  certain  extent  can  be  identified  with  surfaces  of  constant  pressure 
[8],  and  in  color  of  the  local  density.  In  the  fourth  column  the  pressure,  density,  temperature,  Mach 
number  Maz ,  and  velocity  w  component  in  the  z-direction  along  the  axis  of  the  vortex  are  shown. 
The  pictures  clearly  indicate  that  with  increasing  time  a  free  stagnation  point  on  or  near  the  axis  is 
formed,  followed  by  a  bubble-like  flow  structure.  The  bubble  grows  in  the  lateral  and  axial  directions 
in  time,  until  a  stable,  only  slightly  oscillating  state  is  reached.  Inside  the  bubble  several  small  ring-like 
vortex  structures  can  be  identified,  traveling  downstream  with  progressing  time.  From  the  velocity 
and  Mach  number  distributions  it  can  be  seen  that  there  is  a  strong  axial  reverse  flow  along  the  axis 
of  the  vortex.  The  A2-surface  representations  confirm  the  ring-like  structures  inside  the  bubble.  For 
increasing  time  the  flow  deviates  more  and  more  from  its  initial  axial  symmetry,  confirming  that  it  can 
be  concluded  that  vortex  breakdown  is  a  truly  three-dimensional  process  which  cannot  be  described 
using  a  formulation  for  axially  symmetric  flows. 

The  schlieren  picture  in  Fig.  7,  shows  the  bubble  topology  of  the  computation.  In  the  parallel  exper¬ 
imental  investigations  the  vortex  is  produced  with  a  double  wedge  airfoil  in  a  supersonic  flow.  When 
vortex  breakdown  occurs,  the  shock  is  deformed  upstream,  and  the  characteristic  length  scale  of  the 
arising  bubble  is  of  the  order  of  the  diameter  of  the  core  of  the  impinging  vortex. 

In  addition  to  the  numerical  and  experimental  investigations,  a  breakdown  criterion  was  derived  for 
uniform  axial  flow  (Fig.  8).  Based  on  the  axial  momentum  equation  for  inviscid  flow  and  the  Rankine- 
Hugoniot  relations,  the  onset  of  breakdown  can  be  predicted  by  requiring  a  stagnation  point  to  be 
formed  on  the  axis.  In  reference  [2]  it  is  shown  that  a  breakdown  condition  based  on  the  momentum 
equation  along  the  vortex  axis  reads 
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In  the  above  equation,  the  pressure  and  momentum  forces  are  written  with  respect  to  the  free  stream 
pressure  p*^ .  In  the  case  of  uniform  axial  flow,  the  second  term  (II)  on  the  left  side  can  be  written 
in  the  following  form 

n*  .  7,2* 

AXIS  vz1axis 


Since  CROCCOS  equation 

TVs  +  fx(Vxf)  =  ^+  V/io  (4) 

yields  T  •  Vs  =  V/io  if  the  velocity  and  the  vorticity  vector  are  parallel  to  each  other  in  steady  flows 
the  radial  momentum  equation  for  slender  vortices 
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can  be  solved  for  the  velocity  distribution  given  by  equation  (1).  Together  with  the  Rankine  Hugoniot 
relation  for  the  pressure  ratio  across  a  normal  shock  (term  III) 


2-7 
7  +  1 


(Mai  ~  1) 
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the  axial  momentum  equation  (2)  leads  to  the  following  breakdown  criterion 


(1 _ _ )* 

167t2[2  +  (7  —  1  )Ma%0]~1 


jMa2lo 


(1- 


(7— i)r^-ei 


167r2[2+(7-l)MaJJ 


—  U"1 


>  1  + 


27  (Mai-  1) 


7+1 


(7) 


Pi  4 


P\A 


Jp*  o 


P2oo/Pl0 


From  the  above  equation  the  circulation  r0  for  which  breakdown  occurs  can  be  determined.  The 
corresponding  spiral  angle  r  is  defined  by  the  quotient  of  the  circumferential  and  the  axial  velocity 


7(r*  =  tq)  =  r°  •  \A  +  +Ma^°  (8) 

T  uj(r*  =  Tq)  2  •  7r  •  Maoo  •  (1  —  S  •  e~^w) 

as  a  function  of  the  free  stream  Mach  number  Ma^.  The  comparison  of  the  results  obtained  using 
the  above  relation  with  theoretical,  numerical  and  experimental  data  [9,  10,  11]  shows  good  agreement 
in  the  range  of  1.5  <  Ma^  <  2.0  (Fig.  9). 


4  Concluding  Remarks 

Breakdown  of  longitudinal  slender  vortices  in  supersonic  flow,  caused  by  the  interaction  with  a  normal 
shock,  was  investigated  using  a  numerical  solutions  of  the  Euler  and  Navier- Stokes  equations.  The 
investigation  of  the  normal  shock- vortex  interaction  yields  two  types  of  flow  fields  downstream  from 
the  shock.  Depending  on  the  inflow  parameters  a  weak  interaction  with  no  vortex  breakdown  was 
observed,  and  a  strong  interaction  with  a  bubble-like  vortex  breakdown.  The  calculations  show,  that 
a  slender  vortex  propagating  across  a  normal  shock  does  not  necessarily  leads  to  breakdown.  If  the 
circulation  of  the  vortex  is  small  and  if  the  axial  velocity  profile  is  uniform,  the  vortex  does  not  change 
its  overall  shape.  By  increasing  the  vortex  strengths,  the  free  stream  Mach  number  or  the  vortex 
deflect  of  the  axial  flow  the  shock  vortex  interaction  leads  to  a  free  stagnation  point  on  or  nearby  the 
vortex  axis,  which  can  be  understood  as  the  initialization  of  vortex  breakdown.  The  results  obtained 
indicate,  that  the  flow  structure  of  the  burst  part  of  the  vortex  as  simulated  with  the  numerical  solution 
agrees  well  with  experimental  observations.  Several  ring-like,  slightly  oscillating  vortex  structures  are 
formed  immediately  downstream  from  a  stagnation  point  on  the  axis.  For  certain  flow  conditions 
strong  upstream  flow  motion,  even  supersonic  flow  may  occur  near  the  axis.  A  breakdown  criterion 
was  derived  that  showed  good  agreement  with  numerical  and  experimental  findings. 
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Figure  1:  Computational  setup  and  boundary  conditions  for  normal  shock- vortex  interaction.  Flow 
is  in  the  positive  ^-direction. 
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Figure  2:  Vortex  parameters  as  a  function  of  r  at  Maoo  =  1.6,  7=1.4,  ro  =  1.0..5.0 
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Figure  3:  Mach  number  distribution  (left  column)  and  streamlines  (right  column)  in  the  midplane 
of  the  computational  domain  at  Maoo  =  1.6,  r0=3.5,  4.0,  4.5,  and  inviscid,  uniform  flow  (time  step 
£=200). 
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Figure  4:  Instantaneous  pictures  at  time  step  t— 200  of  the  Euler  (left  column)  and  Navier-Stokes 
(right  column)  computations  (Maoo  =  1.6,  ro=5.0,  and  (5=0).  a)  spatial  iso-vorticity  surface  combined 
with  the  Mach  number  distribution  projected  on  the  lower  boundary,  b)  streamline  pattern  in  the 
x-z- middle  plane  of  the  bubble,  and  c)  Mach  number  distribution  within  the  bubble  at  the  x-y- plane 
for  z  =  0. 
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Figure  5:  Evolution  in  time  of  streamlines  (1st  column),  and  vorticity  contours  (2nd  column)  in  the 
center  plane  y  —  0;  A2-iso  surface  (3d  column);  distribution  of  the  primitive  variables  and  the  Mach 
number  on  the  centerline  (4th  column  from  left  to  right)  for  certain  time  steps  for  normal  shock- vortex 
interaction  (Ma^  =  1.6,  r0  =  2.5,  S  =  0.1). 
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Figure  6:  Flow  schematic  for  strong  vortex  breakdown  at  time  step  t  =  200. 


Figure  7:  Color  schlieren  picture  (M.  Klaas,  private  communications)  (left)  and  numerical  Mach 
number  distribution  (right)  at  Ma^  =  1.6,  r0=4.5,  (5=0  for  uniform  axial  flow  at  i?e=15000.  Flow  is 
from  left  to  right. 


Figure  8:  Theoretical  model  to  predict  vortex  breakdown. 
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Figure  9:  Breakdown  map  for  normal-shock  vortex-interaction. 
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Question  by  Mr.  Jeune :  Concerning  the  vortex  breakdown,  limit  presented  as  a  function  of  Mach 
number  for  normal  shock  case,  should  it  be  possible  to  correlate  this  limit  for  both  normal  and 
oblique  shocks  by  considering  the  pressure  jump  across  the  shock? 

Answer:  No,  unfortunately  not.  In  case  of  normal  shocks  the  flow  downstream  of  the  pressure 
jump  is  subsonic.  This  is  a  “must”  for  the  decelerated  flow  on  the  vortex  axis  in  order  to  obtain 
a  free  stagnation  point.  In  case  of  an  oblique  shock,  not  only  must  the  pressure  jump  across  the 
shock  be  a  criterion,  but  also  the  “strong”  shock  solution  must  be  obtained  locally  in  the  vicinity 
of  the  vortex  core.  Therefore,  locally  at  vortex  axis,  a  small  subsonic  pocket  has  to  be  formed 
in  which  the  flow  is  decelerated  so  that  a  stagnation  point  becomes  possible.  By  taking  a  look 
on  the  local  Mach  numbers,  the  water  parameters  5  can  be  varied  to  reduce  the  Maxis  such  that 
Maxis  «  M00.  The  flow  on  the  vortex  axis  is  not  able  anymore  to  follow  the  redirection  (given 
by  the  ramp  angle)  so  that  locally  a  normal  shock  is  formed.  If  now,  as  a  second  criterion,  the 
pressure  jump  across  this  locally  normal  shock  is  sufficient,  the  vortex  will  break  down. 

For  further  information  about  that,  please  take  a  look  at  my  ECCOMAS  paper  2000,  listed  in  the 
references. 
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